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Abstract 



The charged-particle pseudorapidity density dNch/drj has been measured 
for Au+Au collisions at ^/s]^=130 GeV at RHIC, using the PHOBOS ap- 
paratus. The total number of charged particles produced for the 3% most 
central Au+Au collisions for |r/| <5.4 is found to be 4200+470. The evolution 
of dNch/df] with centrality is discussed, and compared to model calculations 
and to data from proton-induced collisions. The data show an enhancement 
in charged-particle production at mid-rapidity, while in the fragmentation re- 
gions, the results are consistent with expectations from pp and pA scattering. 
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Nucleus-nucleus collisions at highly relativistic energies from the Relativistic Heavy Ion 
Collider (RHIC) present a new opportunity for the study of strongly interacting matter 
under conditions of very high temperature and energy density. Already, data from central 
Au+Au collisions at energies of ^/sNN=5Q and 130 GeV have shown an increase in 
mid-rapidity particle production over that observed in Pb+Pb collisions at the CERN SPS 
0. In addition, the charged-particle yields, when scaled by the number of participating 
nucleons (Npart), exceed the values observed in proton-proton collisions at similar s/snn- 
These results have been used to infer an energy density at the highest value of ^/snn that is 
considerably larger than was achieved at CERN, and which is well within the region where 
new physics is to be expected. 

The distribution of dNch/di] over the entire range of pseudorapidity rj, where i] = 
—ln{ta.n{6/2)) and 6 is the angle with respect to the beam, represents a time-integral of 
particle production throughout the collision, and reflects all effects that contribute to the 
production of charged particles. In addition to the influence of both hard and soft produc- 
tion processes, dNch/dr] is sensitive to nuclear effects in the initial parton distributions, as 
well as to effects from hadronic re-interactions in the final state. 

The data for dNch/drj at mid-rapidity (|?7| < 1) for the most central Au+Au collisions at 
\/snn=56 and 130 GeV |]1| are in reasonable agreement with the predictions of a number of 
models including HIJING 0, a saturation model(EKRT) [Q, and purely hadronic models 
(e.g. LUCIFER PJ^). The centrality dependence of dNch/dr] at mid-rapidity has also been 
measured These results suggest some deviation from both the HIJING and EKRT 

model predictions, although they are in broad agreement with the results of calculations by 
Kharzeev and Nardi [^]. 

An extension of the measurements of dNch/dr] data beyond mid-rapidity, for a range of 
impact parameters, is necessary to further constrain models. It is also of interest to determine 
whether the observed scaling of the charged-particle yield with N^art at mid-rapidity is 
modified at large values of r] where, in proton-nucleus {pA) collisions [pUj-P^, re-scattering, 
stopping, and target fragmentation influence the shape of the dN^h/dr] distributions. In 
order to address these questions, we have used the PHOBOS apparatus to measure the 
charged-particle pseudorapidity density dNch/dr] from Au+Au collisions at i/s5v7v=130 GeV 
over the range \r]\ < 5.4. 

The PHOBOS experiment at RHIC largely consists of several arrays of silicon pad detec- 



tors. The details of the experimental arrangement are described in Ref. [|13|. The procedures 



used for event selection, the determination of the collision-vertex position, event centrality, 
and the estimation of Npart, have been described in Refs. The specific elements of the 

experiment used in the current measurement and the analysis procedures leading to dNch/dr] 
are described below. 

The data samples included in the current analysis were taken at a collision energy of 
y^sjv7v=130 GeV. The collision vertices were confined to a region within ±10 cm of the 
nominal beam crossing and center of the apparatus {z = 0). At mid-rapidity, charged 
particles were detected, and their energy deposition measured, with an octagonal array 
of pad detectors approximately 1 m long (the "Octagon") that surrounds the thin-walled 
Be beam pipe. The Octagon subtends the full azimuthal range, except for regions where 
sensors that would intercept the acceptance for the tracking spectrometers and vertex finding 
detectors are removed. For collision vertices within \z\ <10 cm, the pseudorapidity coverage 
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of the Octagon is complete for |?7| <3.2. Six rings of silicon pad detectors placed at distances 
of \z\ ~1, 2 and 5 m, detected particles in pseudorapidity ranges of 3 < \'q\ < 4, 4 < \rj\ < 4.7, 
and 4.7 < |?7| < 5.4, respectively. Thus, for collisions within the 20 cm long region in the 
center of the experiment, there are no significant gaps in the rj coverage from one sub- 
detector to the next. The total numbers of pads in the Octagon and Rings are 11,040 and 
3,072, respectively. 

Two complementary methods were used to analyze the pseudorapidity distribution data. 
The first method ("Hit Counting") uses the segmentation of the multiplicity detector. After 
merging of signals in neighboring pads, in cases where a particle travels through more than 
a single pad, the deposited energy was corrected for the angle of incidence, so that all 
tracks originating from the collision vertex possess a common average value of the deposited 
energy {AE k, 80 keV). Pads containing more than 75% of this value were counted as 
occupied. This requirement largely suppresses hits from background and from secondary- 
particles not originating from the primary collision vertex. Then, for a given value of r] and 
bin i in collision centrality, the number of hit pads N{ri, i) was corrected for the effects of 
multiple occupancy, where more than one particle travels through a given pad, as well as 
for contributions from the remaining secondary particles, absorption in the beam pipe, and 
weak decays of primary particles. The number of particles per pad was determined as a 
function of rj and centrality bin in two ways. In the first method, the probability P{N) of 

particles passing through a given pad was assumed to be Poisson distributed. The mean 
occupancy /i(?7, i) could then be determined from the ratio of occupied to unoccupied pads in 
a range of i] for each centrality bin. Second, as a check, the occupancy was also determined 
from the energy-deposition spectra. A fitting procedure was used to determine the relative 
contribution to the total energy deposition of one or more particles, thus providing an 
independent measure of the mean occupancy ^{ri,i). Maximum values of approximately 1.6 
particles per hit pad were obtained for the most central collisions at mid-rapidity. 

To account for effects or biases not treated in the above analysis procedures, a final 
correction was deduced from GEANT simulations of the detector response using events 
from the HIJING g, RQMD H] and VENUS event generators. The ratio between the 
simulated, occupancy-corrected dNch/drj distributions and the known "Truth" distributions 
formed this last set of corrections, which ranged up to 15% in the Octagon and up to 50% in 
the Rings. These background correction factors, dependent upon both rj and centrality bin, 
were applied to the occupancy-corrected data, yielding dNch/drj. The final results obtained 
using both hit counting methods were in good agreement. 

In the second ( "Analog" ) method, the pseudorapidity distribution was extracted directly 
from the energy deposition AE{r]) in the multiplicity detectors. The measured energy depo- 
sition was transformed into dNch/dr) using quantities derived from the results of Monte-Carlo 
simulations. The average energy per track {AEtr) was determined as a function of rj using 
particles from HIJING events, passed through the GEANT simulation of the detector. The 
fraction of primary particles fprimiv) determined from the same simulations. Then, 
dNch/dr] = ^^^E^ ^ri)Arj'^^ " Although the two methods differ qualitatively in the way in which 
the energy-deposition information from the multiplicity detectors is used, and rely differ- 
ently upon the results of Monte-Carlo simulations, they yield results that differ generally by 
^5% throughout the range in t], well within the systematic uncertainties (see below), which 
are approximately 10%. 
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The systematic uncertainty in the occupancy correction for the hit-counting analysis 
was obtained by comparing the results from the full analysis chain using Poisson-derived 
occupancy corrections with those derived from the measured AE spectra. The average 
deviations are less than 3%, yielding a partial systematic error of ~ 3% at mid-rapidity. 
The systematic uncertainties from the Monte-Carlo simulations have been estimated by 
using different assumptions in the GEANT simulation, as well as different event generators, 
including RQMD and VENUS The variations observed in the derived background 



corrections are between 4 and 8%, suggesting a total systematic uncertainty of approximately 
10%. The systematic errors in the Analog analysis arise from uncertainties in the Monte 
Carlo simulations, as well as those in the absolute energy calibrations of the Silicon-pad 
detectors. The latter are approximately 5%, yielding a total systematic error for the Analog 
method of approximately 10%, similar to that for the hit-counting method. 

Our final results are presented in Fig. |l|(a)-(f), which show the error-weighted average 
values of dNch/dt] from the two procedures for six different centrality bins. The error bars 
represent a convolution of the estimated systematic errors in the different analyses. The 
different centrality bins are denoted by the corresponding fraction of the observed total cross 
section, as well as by the deduced average number of participant nucleons < Npart >■ While 
the occupancy and background corrections applied to the data from different parts of the 
apparatus are quite different in magnitude, the data from the Octagon and Ring detectors 
merge smoothly in most cases. The present analyses give values of dNch/drj that are in good 
agreement from the independent "Tracklet" analysis presented in P). For example, for the 
6% most central collisions, we find dNch / dri\\n\<i = 547±55, compared to 580±25 from Ref. 

i- 

The integral of the distributions N^^, plotted in Fig. ||(a) as a function of centrality, is a 
direct measure of the total entropy produced in the collisions. Predictions of its magnitude 
have varied by as much as factor of two fl^. With increasing (Npart), the observed values 
of iV*^* change smoothly from 910 ± 100 for the 40-45% centrality bin to 4200 ± 470 for the 
3% most central collisions. Per participant pair, these numbers correspond respectively to 
21.8±2.6 and 23.7±2.7, compared with the pp and pp non-diffractive total charged particle 
multiplicity of 18.5±0.7 [0. The predictions of the HIJING model reproduce the general 



trend of the A^*)J* centrality dependence, but systematically under-predict the observed values 
by approximately 10%. 

The shapes of the pseudorapidity distributions evolve gradually with increasing cen- 
trality, as shown in Fig. |I|(a-f). For all centrality bins, there is a plateau region between 
— 2 < < 2, followed by a rapid drop-off towards larger pseudorapidities. A more detailed 
study of the centrality dependence of the shape is given in Figs. 2(b-f) and Fig 3(a). In 
Figs. 2(b-f), the Npart dependence of dNch/drj normalized per participant pair < Npart /'^ >, 
is plotted for five pseudorapidity bins ranging from |?7| < 1 to 5 < \7]\ < 5.4. Also plotted 
are data from pp and pp [|19[ collisions, scaled as described below, as open circles, and 



predictions from the HIJING model as solid lines. The statistical uncertainties are small and 
the systematic uncertainties are comparable to those described above. Figure 3(a) shows the 
pseudorapidity distribution for peripheral (35-45%) and central (0-6%) Au+Au collisions, 
scaled by the respective number of participants. 

For all pseudorapidity bins in Fig. 2(b-f), the data evolve smoothly from the most pe- 
ripheral to the most central collisions. As seen in Fig. 2(b,c) and Fig. 3(a), central collisions 
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yield a 10-15% higher dNch/di] per participant in the plateau region, compared to periph- 
eral events. This difference decreases at larger pseudorapidities and reaches a turnover point 
between |?7| = 3 and 4, as seen in Fig. 3(a). Beyond that, a higher yield per participant is 
found in peripheral collisions. In Fig. 2(d) and (e) this turnover can be seen as a change in 
slope of the Npart dependence of dN/dt] in the respective t] bins. In the highest 1.5 units 



of pseudorapidity, as is observed in p + A collisions at lower energies iToHr^, the scaled 
charged-particle density actually falls with Npart, being reduced by nearly a factor of two for 
|?7| > 5. Qualitatively, the changes in the distributions from peripheral, to central collisions 
in the fragmentation regions are similar to those observed in pA collisions |]TD|-0 . 

Finally, a comparison of the shapes from pp/pp collisions and central Au+Au collisions 
is given in Fig. 3(b). It shows the scaled pseudorapidity distribution for central collisions 
(0-6%), compared to scaled data from pp and pp collisions, shown as a grey band. The 
pp/pp distribution was obtained by scaling the measured dN/drj distributions from pp col- 
lisions at -y/i = 53 GeV and pp collisions at ^/s = 200 and 546 GeV [l^Jl^ horizontally by 



2/max(130 GeV) / ymaxiV^) , whcrc VmaxiVs) = ^^iVs/nip), and vertically using the parame- 
terization of dN/dr]\rj\^i from [|1^[|. The grey band gives an estimate of the uncertainty of 
this extrapolation procedure. 

The scaled dNch/dr] is found to be higher in central Au-|-Au collisions than in pp/pp over 
the full pseudorapidity range out to |?7| > 4, with the largest excess observed in the central 
plateau region. This is in contrast with the HIJING prediction, which shows an excess in 
Au-|-Au collisions only for \ri\ < 2 — 3. A possible origin of this qualitative difference is 
suggested by the AMPT model of Zhang et al. |2^. In this model, the initial state parton 
distribution is obtained the same way as in HIJING, but is followed by a parton cascade 



|21|, string fragmentation and hadronic re-scattering using a relativistic transport model 
2^. This model reproduces the excess in particle production at higher pseudorapidities 
seen in Au-|-Au collisions relative to pp/pp. Furthermore, predictions of dNch/drj in a com- 
pletely hadronic framework (e.g. LUCIFER) are very similar to those of the AMPT 
calculations, and are also in good agreement with the data for the most central collisions. 
These observations suggest that effects in the hadronic phase should be taken into account to 
provide a full description of the data. A comparison of the predicted centrality dependence 
of the full distributions from these models with the data would also be of interest. Further 
insight will be gained from future RHIC data, which should include more precise reference 
data from pp collisions, as well as nucleus-nucleus data at different collision energies. These 
results will provide a basis for separating the effects of new phenomena from conventional 
hadronic physics. 

This work was supported in part by US DoE grants DE-AC02-98CH10886, DE-FG02- 
93ER-404802, DE-FC02-94ER40818, DE-FG02-94ER40865, DE-FG02-99ER41099, W-31- 
109-ENG-38, and NSF grants 9603486, 9722606, and 0072204. The Polish groups were 
partially supported by KBN grant 2P03B 04916. The NCU group was partially supported 
by the NSC of Taiwan under contract NSC 89-21 12-M-008-024. 



^This empirical scaling is chosen such that the extrapolation of three sets of pp or pp data are 
confined within the band shown on Fig. ^(b). 
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FIG. 1. Charged-particle pseudorapidity density dNch/drj from y^sjvAf=130 GeV Au+Au col- 
lisions, for different centrality bins, as defined by different fractions of the total observed cross 
section. The error bars reflect primarily the systematic uncertainties. The average number of 
participants < Npart > for each bin is also indicated. For fractions of the cross section > 45% 
the systematic uncertainties in the < Npart > determination are still under study and no value is 
quoted. 
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FIG. 2. (a): Total number of charged particles detected within the range —5.4 < rj < 5.4. 
The symbols are the measured data, and the line represents the prediction of the code HIJING. 
The error bars reflect systematic uncertainties, (b)-(f): Centrality dependence of dNch/drj for 
different ranges of |?7|. The filled symbols are the measured data, and the solid curves are the 
HIJING predictions. The open diamonds and triangles refer to the predictions of the AMPT, and 
LUCIFER, models, respectively for the 6% most central collisions. The open circles represent the 
values from pp and pp collisions. 
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FIG. 3. (a). dNch/dr]/{< Npart > /2) for central (0-6%, open symbols) and peripheral (35-45%, 
filled symbols) collisions. The error bars have been suppressed for clarity, (b) Comparison of the 
measured pseudorapidity distribution for the 0-6% centrality bin with the HIJING (thin curve) 
and AMPT (dashed curve) models. The grey band represents data from pp and pp collisions 
interpolated to y/sNN=l30 GeV. 
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